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Abstract 



Theoretical expectations concerning small x behaviour of the spin dependent 
[ structure function gi are summarised. This includes discussion of the Regge pole 

model predictions and of the double ln 2 (l/x) effects implied by perturbative 
I> \ QCD. The quantitative implementation of the latter is described within the uni- 

fied scheme incorporating both Altarelli-Parisi evolution and the double ln 2 (l/x) 
resummation. The double ln 2 (l/x) effects are found to be important in the re- 
gion of x which can possibly be probed at HERA. Predictions for the polarized 
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gluon distribution AG(x,Q ) at low x are also given. 



Understanding of the small x behaviour of the spin dependent structure functions of 
the nucleon, where x is the Bjorken parameter is interesting both theoretically and 
phenomenologically. Present experimental measurements do not cover the low x val- 
ues and so the knowledge of reliable extrapolation of the structure functions into this 
region is important for estimate of integrals which appear in the Bjorken and Ellis- 
Jaffe sum rules pj. Theoretical description of the structure function gi(x,Q 2 ) at low 
x is also extremely relevant for the possible polarised HERA measurements [0. The 
purpose of this talk is to summarize the theoretical QCD expectations concerning the 
small x behaviour of the spin dependent structure function gi of the nucleon. After 
brief reminder of the Regge pole expectations we shall discuss the effects of the double 
ln 2 {l/x) resummation and its phenomenological implementation. Besides the structure 
function g x we shall also discuss the spin dependent gluon distribution AG. 
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The small x behaviour of spin dependent structure functions for fixed Q 2 reflects 
the high energy behaviour of the virtual Compton scattering (spin dependent) total 
cross-section with increasing total CM energy squared W 2 since W 2 = Q 2 (l/x — 1). 
This is, by definition, the Regge limit and so the Regge pole exchange picture [|J is 
therefore quite appropriate for the theoretical description of this behaviour. Here as 
usual Q 2 = —q 2 where q is the four momentum transfer between the scattered leptons. 
The relevant Reggeons which describe the small x behaviour of the spin dependent 
structure functions are those which correspond to the axial vector mesons 0. [J| . 

The Regge pole model gives the following small x behaviour of the structure func- 
tions g\(x, Q 2 ) where g\(x, Q 2 ), i = s,ns denote either singlet (gf(x, Q 2 ) = g\(x, Q 2 ) + 
g™ (x, Q 2 )) or non-singlet (g™ s (x, Q 2 ) = gf(x, Q 2 ) — gi(x, Q 2 )) combination of structure 
functions: 

g\{x,Q 2 )= ll {Q 2 )x- a ^ (1) 

where g\(x, Q 2 ) denote either singlet {g({x, Q 2 ) = gi(x, Q 2 ) + g™(x, Q 2 )) or non-singlet 
(gi"(x, Q 2 ) = gl(x, Q 2 ) — g™(x, Q 2 )) combination of structure functions. 

In eq. and flU) ot s>ns {Q) denote the intercept of the Regge pole trajectory correspond- 
ing to the axial vector mesons with I = or I = 1 respectively. It is expected that 
a s ,ns(0) < and that a s (0) ~ a ns (0) i.e. the singlet spin dependent structure function 
is expected to have similar low x behaviour as the non-singlet one. 

It may be instructive to confront this behaviour with the Regge pole expectations 
for the spin independent structure function Fi(x, Q 2 ): 

Fl(x,Q 2 )=(3 l (Q 2 )x- a ^. (2) 

The singlet part Ff = Ff + F™ of the structure function F\ is controlled at small x 
by pomeron exchange, while the non-singlet part F^ s = Ff — F™ by A2 reggeon. 
The pomeron intercept is significantly different from that of the A2 Reggeon i.e. 
ap(0) = l + e with e ~ 0.08 while 0:^2(0) ~ 0.5 as determined from the fits to the total 
hadronic and photoproduction data || . This implies that in the spin independent case 
the singlet part Ff(x, Q 2 ) of the structure function F\(x, Q 2 ) dominates at low x over 
the non-singlet component. 

Several of the Regge pole model expectations for both spin dependent and spin in- 
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dependent structure functions are modified by perturbative QCD effects. In particular 
as far as the spin dependent structure functions are concernced the Regge behaviour (II]) 
becomes unstable against the QCD evolution which generates more singular behaviour 
than that given by eq. ([!]) for «i(0) < 0. In the LO approximation one gets: 

g? s (x,Q 2 ) ~ exp[2^AP qq m(Q 2 )ln(l/x)} 



g*(x, Q 2 ) ~ exp[2^n(Q 2 )Hl/x)} (3) 

where 

f«?) = ^f (4) 
JqI q 2 2n 

and 



AP qq (0) + AP ag (0) + ^(AP qq (0) - AP 99 (0)) 2 + 4AP (?9 (0)AP OT (0) 
7 + = g 

with APjj(O) = APij(z = 0) where APy (z) denote the LO splitting functions describing 
evolution of spin dependent parton densities. To be precise we have: 

AP 99 (0) = ^ (6) 
AP w (0) = -N F 

AP M (0) = | 
AP 33 (0) = 12. 

We recall for comparison the analogous small x behaviour of the structure function 
Fi(x, Q 2 ) in the LO approximation : 



Pf(x,Q 2 ) ~ -exp[2^/6aQ 2 )ln(l/x)}. (7) 

The Regge behaviour of the non-singlet structure function F± (x,Q 2 ) remains stable 
against the QCD evolution. 



The LO (and NLO) QCD evolution which sums the leading (and next-to-leading) 
powers of IuQ 2 /Qq) is however incomplete at low x. In this region one should worry 
about another "large" logarithm which is ln(l/x) and resum its leading powers. In 
the spin independent case this is provided by the Balitzkij, Fadin, Kuraev, Lipatov 
(BFKL) equation which gives in the leading ln(l/x) approximation the following 
small x behaviour of Ff(x, Q 2 ) 

F?( X ,Q 2 ) ~ x~ Xbfkl (8) 
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where the intercept of the BFKL pomeron Xbfkl is given in the leading order by the 
following formula: 

Xbfkl = l + —*ln(2). (9) 

7T 

It has recently been pointed out that the spin dependent structure function g x at 
low x is dominated by the double logarithmic ln 2 (l/x) contributions i.e. by those 
terms of the perturbative expansion which correspond to the powers of ln 2 {l/x) at 
each order of the expansion [H, H. Those contributions go beyond the LO or NLO 



order QCD evolution of polarised parton densities QTOj and in order to take them into 
account one has to include the resummed double ln 2 (l/x) terms in the coefficient and 
splitting functions [n||. In this talk we will present discussion of the double ln 2 (l/x) 
resuumation following alternative approach based on unitegrated distributions |L2L |13| . 



The dominant contribution to the double ln 2 (l/x) resummation comes from the 
ladder diagrams with quark and gluon exchanges along the ladder (cf. Fig.l). In what 
follows we shall neglect for simplicity possible non-ladder bremsstrahlung terms which 
are relatively unimportant || || . 



It is convenient to introduce the unintegrated (spin dependent) parton distributions 
fi(x',k 2 ) (i = u v ,d v ,u,d,s, g) where k 2 is the transverse momentum squared of the 
parton i and x' the longitudinal momentum fraction of the parent nucleon carried by 
a parton. The conventional (integrated) distributions Api(x, Q 2 ) are related in the 
following way to the unintegrated distributions fi(x', k 2 ): 

Ap l (x,Q 2 ) = Ap ( ?\x) + ™ fi ( x ' = x (i+K), k 2 ) (10) 

Jkl k 2 Q 2 

where Apf^ (x) is the nonperturbative part of the of the distributions. The parameter 
&o is the infrared cut-off which will be set equal to 1 GeV 2 . The origin of the nonper- 
turbative part Apf*\x) can be viewed upon as originating from the non-perturbative 
region k 2 < k" 2 , i.e. 

Apf\x)= J^-^f^k 2 ) (11) 

The spin dependent structure function gf (x, Q 2 ) of the proton is related in a standard 
way to the (integrated) parton distributions: 



9 P i(x,Q 2 ) = 

~{Au v {x, Q 2 ) + 2Au{x, Q 2 )) + -(Ad v (x, Q 2 ) + 2Au(x, Q 2 ) + 2As{x, Q 2 )) ( 1 2 ) 
y y 
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Figure 1: An example of a ladder diagram generating double logarithmic terms in the 
spin structure function g 1 . 

where Au v (x,Q 2 ) = Ap Uv (x,Q 2 ) etc. We assume Aw = A J and confine ourselves to 
the number of flavours N F equal to three. 

The valence quarks distributions and asymmetric part of the sea 

f us (x',k 2 ) = f a (x',k 2 )-f- s (x',k 2 ) (13) 

will correspond to ladder diagrams with quark exchange along the ladder. The singlet 
distributions 

fs(x\ k 2 ) = f Uv (x', k 2 ) + f dv (x', k 2 ) + Af a (x', k 2 ) + 2f,(x', k 2 ) (14) 

and the gluon distributions f g (x',k 2 ) will correspond to ladder diagrams with both 
quark (antiquark) and gluon exchanges along the ladder. 

The sum of double logarithmic ln 2 {l/x) terms corresponding to ladder diagrams is 
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generated by the following integral equations (see Fig.l): 

*dz r k2 / z dk' 2 



f k (x\e) = fj?\x>,e) + ^AP qq (o) f ^ / 

Z7T Jx' Z Jk 

(yk —— dy^ us) 



x 



k» fh [7> k 



12 



(15) 



fs(x',k 2 ) = ff\x',k 2 ) + 



a s r^dz r k2 /*dk' 2 
2tt Jx' z Jkl k' 2 



AP qq (0)fslj,k' 2 \+AP qg (0)f g i^k 



12 



f g (x',k 2 )=6°\x',k 2 ) + 



l dz r k2 / z dk 



12 



AP gq (0)f s i^k' 2 \ + AP gg (0)f g i^,k' 2 



(16) 



27T Jx' Z Jk^ k' 2 

with APij(O) = APij(z = 0) given by eq. (0). 

The variables k 2 (k' 2 ) denote the transverse momenta squared of the quarks (gluons) 
exchanged along the ladder, k 2 is the infrared cut-off and the inhomogeneous terms 
fi (x', k 2 ) will be specified later. The integration limit k 2 / z follows from the require- 
ment that the virtuality of the quarks (gluons) exchanged along the ladder is controlled 
by the tranverse momentum squared. 



The origin of the double logarithmic ln 2 (l/a;) terms in gi(x, Q 2 ) can be traced to 
the fact that the conventional single logarithmic terms coming from the logarithmic 
integration over the longitudinal momentum fraction z are enhanced by the logarith- 
mic integration over the transverse momentum up to the z dependent limit k 2 / z in 
equations (|I5||T6D and up to the x dependent limit W 2 = Q 2 {l/x — 1) in eq. (ITT 



Equation flT5| ) is similar to the corresponding equation in QED describing the dou- 
ble logarithmic resummation generated by ladder diagrams with fermion exchange WA 



The problem of double logarithmic asymptotics in QCD in the non-singlet channels was 



also discussed in ref 15 



Equations (15], [16]) generate singular power behaviour of the spin dependent parton 
distributions and of the spin dependent structure functions g± at small x i.e. 



9i~(x,Q ) ~ x 
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AG(x, Q 2 ) ~ x~ Xs (17) 

where g± s = gf — g™ and gf = g^ + g™ respectively and AG is the spin dependent gluon 
distribution. This behaviour reflects similar small x' behaviour of the unintegrated 
distributions. Exponents Xns,s are given by the following formulas: 



A 



NS 



Ac 



2^ 



AP M (0) 



2^ 



7" 



where 7 + is given by eq.(|5]). The power-like behaviour ([17]) with the exponents Xns,s 
given by eq. (|18D remains the leading small x behaviour of the structure functions pro- 
vided that their non-perturbative parts are less singular. This takes place if the latter 
are assumed to have the Regge pole like behaviour with the corresponding intercept (s) 
being near 0. 



In order to understand origin of the power-like behaviour ( |TTD it is useful to go to 
the moment space and inspect the singularities of the moment functions fi(w, k 2 ) and 

APi(LU,Q 2 ) 



fi(u;,k 2 )= / dx'x'^f^x'X) 
Jo 

Api(u,Q 2 )= f 1 dxx^Ap^x^ 2 
Jo 



(19) 

in the u plane. It follows from eq. ([T0|) that the moment functions Api(u>,Q 2 ) are 
related in the following way to fi(u>, k 2 ): 



a Pj ( W) q 2 ) = a#h + 



00 dk 2 ( k 2 
1 + 



kl k 2 



Q 2 



fi(u,k 2 



(20) 



where Ap\°\u) denote the moment functions of the nonperturbative distributions 
Apf\x). Let us first consider the d v ,us. Equation (|15|) implies the 

following equation for the moment functions fi(w, k 2 ) 



f t (oj,k 2 ) = f?\u;,k 2 ) + 



In this equation ft (u, k 2 ) denote the moment functions of f\ u \x\ k 2 ) and a s is defined 
by: 

a s = AP qq (0)-^- (22) 

Z7T 



O.s 



CO 



k 2 

M^,k' 2 ) + 



l k' 2 



dk' 2 ( k 2 



k 2 k' 2 \ k 



r> 



>2\ 



(21) 
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Equation (j2TD follows from (p~5|) after taking into account the following relation: 

(23) 



o z \k' 2 I io 



Q(k 2 -k' 2 ) + (^\ Q(k' 2 -k 2 ) 

For fixed coupling a s equation can be solved analytically. Assuming for simplicity 
that the inhomogeneous term is independent of k 2 (i.e. that f^°'(oj,k 2 ) = Ci(u) ) we 
get the following solution of eq.(^TJ): 

/^,£; 2 ) = ^(a^)(^J (24) 

where 

^ NS {a s ,uj) = (25) 

and 

Ri(a s ,u) = CAou) . (26) 

a s 

Equation (|25"1) defines the non-singlet anomalous dimension in which the double log- 
arithmic ln 2 (l/x) terms i.e. the powers of have been resummed to all orders. It 
can be seen from (p5]) that this anomalous dimension has a (square root) branch point 
singularity at oj = Xns 

X NS = 2v^I. (27) 

This singularity will of course be also present in the moment functions fi(u,k 2 ) and 
Api(uj,Q 2 ). It generates singular power - like behaviour of the non-singlet structure 
function gf s (x,Q 2 ) (cf. eq. flTTj) ) . For Ci(u) = — Apf^^uj) the moment functions 
Api(u, Q 2 ) can be shown to have a familiar RG form 

Ap l (uj 1 Q 2 ) = R l (u,a s )r^\ + 0(-^) (28) 

where 

i2i(w,a a ) = Ap> ; (o;) — (29) 

r(w) 

It may be seen from eq. ( ^9|) that the singularity at uo = Xns which is present in the 
anomalous dimension 7^5 (a s , 10) does also appear in the functions Ri . It is therefore 
present in the moment functions Api(u, Q 2 ) for arbitrary values of the scale Q 2 includ- 
ing Q 2 = k 2 . The situation in this case is similar to that in the solution of the BFKL 
equation ]T5|. 



The formula for the exponent A5 which controls the small x behaviour of gf (x, Q 2 ) 
and of AG(x, Q 2 ) can be obtained by inspecting the singularities in the uj plane of 
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the moment functions fs{w, k 2 ) and J g {oJ, k 2 ) which are generated by the correspond- 
ing equations for those functions. They can be obtained from eqs.fllfip and have the 
following form: 



f s (u,k 2 ) = f s \u;,k 2 ) + 



a. 



00 dk 



12 



2-kuj Jkl k' 2 



Q(k 2 - k' 2 ) + Q(k' 2 - k 



k' 2 



AP qq (0)f s (uj, k' 2 ) + AP qg (0)f g (u, k 



/2^ 



/ J (^ 2 ) = /iVi 2 ) + 



a, r°° dk 



n 



2-kuj Jkl k' 2 



Q(k 2 -k' 2 ) + Q(k ,2 -k 2 ) 



k' 2 



AP gq (Q)f s (u, k ) + AP gg (0)f g (uj, k 



/2\ 



(30) 



For k 2 independent inhomogeneous terms (fflioo^k 2 ) = Cg %g (u})) these equations 



have the following solution: 



fs, g (u;,k 2 ) = Rt(u,a s ) 



'k_ 

Ml, 



2\ l + {u,a s ) 



+ R S (u,a 3 ) 



k 

Mo 



2 \ 7 (<•>,£*«) 



(31) 



The (singlet) anomalous dimensions 7 (uj,a s ) are given by the following equations: 



7 (w,a fl ) 



where 7 + is defined by equation (|5]) while 7 is: 



(32) 



AP gg (0) + AP M (0) - ^(AP qq (0) - AP gg (0)) 2 + AAP qg (0)AP gq (0) 
7 = " g ^ 33 ^ 

The functions Kg (u, a s ) can be expressed in terms of Cs j9 {u) and the anomalous 
dimensions 7 (o;,a; s ). The moment functions Aps, g {oJ,Q 2 ) can be shown to have the 
RG form: 



Aps, g (w,Q 2 ) = RjjJu,a a ) f — I + i2£„(w,ar/ ' 



^ g y^,^ SJ , 2 



+ (34) 



It may be seen from eq.(p2[) that the anomalous dimensions 7 ± ( | ^S«s) have (branch 
point) singularities at uj^ 1 

-± 



c^ = 2,/^ 7 ± 



2tt 



(35) 
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with uj + > uj i.e. u> + = As (cf. eqs. flTT|) and (p~8|) ) is the leading singularity. We also 
have X s >> X NS since 7+ >> AP gq (0) (7+ « 1.18 and AP qq (0) = 4/3). This means 
that singlet distributions and singlet structure functions dominate over the non-singlet 
ones at low x. Since this singularity is also present in the functions Rg (u,a s ) it ap- 
pears in Aps t g(uj,Q 2 ) for arbitrary value of the scale Q 2 . It should be noticed that 
the effect of the double ln 2 (l/x) resummation can be quite strong and, in particular 
the exponent As can easily become of the order of unity. This will be seen more ex- 
plicitely in the quantitative implementation of the double ln 2 (l/x) contributions which 
we are going to discuss below. In the exact leading double ln 2 {l/x) approximation the 
anomalous dimensions j NS and 7 ± and the exponents X^s and As acquire additional 
contributions due to bremsstrahlung diagrams. Their effect on Xns,s was estimated in 
refs. II |9| where it was found that the bremmstrahlung terms can enhance Xns by 
about 4% and reduce As by about 10%. 



In order to make the quantitative predictions one has to constrain the structure 
functions by the existing data at large and moderately small values of x. For such 
values of x however the equations fll5|) and (|T6|) are inaccurate. In this region one 
should use the conventional Altarelli-Parisi equations with complete splitting functions 
APij(z) and not restrict oneself to the effect generated only by their z — > part. 
Following refs. [0, [13| we do therefore extend equations (p!5| , |l6D and add to their right 



hand side the contributions coming from the remaining parts of the splitting functions 
APij(z). We also allow the coupling a s to run setting k 2 as the relevant scale. In 
this way we obtain unified system of equations which contain both the complete LO 
Altarelli-Parisi evolution and the double logarithmic ln 2 (l/x) effects at low x. The 
corresponding system of equations reads: 

a s {k 2 )A [idz f k2 /*dk' 2 
3 



fk(x',k 2 ) 
a s (k 2 ) 



[X 



2tt 



fk 

^ dk' 2 4 



+ 



x 

2tt 3 Jx> z Ag k' 2 ~ h {~ 
±dz(z + z 2 )f k (^k' 2 )-2zf k (x',k' 2 ) 



k' 2 \ + 



k' 2 3Jx> z 
a s (k 2 ) fk 2 dk' 2 



2tt 



k' 2 



2 + -ln(l 



1-z 

f k {x\k' 2 ) 



X 



(36) 



f s (x',k 2 ) = f s °\x',k 2 ) + 



as{k 2 
2ir 



±dz r k2 / z dk' 2 A fx 



Z Jk 



12 



a s (k 2 ) fk 2 dk ,2 4 f±dz(z + z 2 )f s (^k' 2 )-2zf s (x',k' 2 ) 



2tt 



k' 2 3 



1-z 



+ 
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a s (k 2 ) [k 2 dk' 2 



a s (k 2 
2ix 



-Nr, 



2tt Jkl k' 2 
1 dz r k2 / z dk 



2+-ln(l-x') 



12 



X 



fs(x',k' 2 ) + 

fc2 dk' 2 r 1 dz 



kl k' 2 . fx' Z 



2zf. 



,x' 



9V _ 5 



f g (x',k 2 ) = fj 3 °\x',k 2 ) + 



(fc 2 ) yidz r*?/*dk 



2tt 



/2 



/s I -,**] + 



z Jkl k' 2 3 J V z 



2tt 



fk 2 dk' 2 r^dz, 4 „ fx' , /2 \ fidz f k2 / z dk' 2 „ fx' , /2 N 



+ 



a s (k 2 ) [ k2 dk' 2 fidz 

' — Qz 



2n 7fc2 k' 2 Jx' z 



f 9 (^k' 2 )-f g (x',_k^_ 2f 



1-z 



+ 



a s (k 2 ) f k2 dk' 2 Ml N F 



2tt Jkl k' 2 



+ 6/n(l — x'y 



fg(X ,k ) . 



(37) 



The inhomogeneous terms ff ) \x\k 2 ) are expressed in terms of the input (inte- 
grated) parton distributions and are the same as in the case of the LO Altarelli Parisi 



evolution 12 



fi°\x',k 2 ) 



a s (k 2 )A ^dz{l + z 2 )Apf{^ 



2zApf(x' 



2tx 3 Jx' z 
ct s (k 2 



2tt 



2 + -ln(l-x') 



1-z 

(0), 



Apf{x' 



(3* 



fP(x',k 2 ) 
a s (k 2 



2tt 3 Jx' z 



a s {k 2 ) 4 fi (1 + z 2 )ApHt) - ZzApgV 



2tt 



1-2 

1 dz 



(2 + -Jn(l - x'))Apf(x') +N F [ —(1 - 2z)Ap<®(- 
3 Jx' 2 y z 



n°\x',k 2 



a s (k 2 ) 
2tt 



I f 1 -(2 - *)Ap<?>(£) + (^ - ^ + 6/n(l - x'))Apf (*' 

O Jx' Z Z Z O 



11 JV* 



a s (£; 2 L A 1 dz 

-D 



27T Jx' Z 



Apf\^)-zApf(x') 



1-z 



+ (1 - 2z)Ap 



(0) 



(39) 



Equations (|36| , |3~7|) together with (^,^) and flT0|) reduce to the LO Altarelli-Parisi 
evolution equations with the starting (integrated) distributions Ap®(x) after we set the 
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Figure 2: Non-singlet part of the proton spin structure function g±(x, Q 2 ) as a function 
of x for Q 2 =10 GeV 2 . Continuous line corresponds to the calculations which contain 
the leading ln 2 (l/x) resummation, , broken line is a leading order Altarelli-Parisi 
prediction, and a dotted one shows the non-perturbative part t g[ NS0,> = #,4/6(1 — x) 3 , 
where qa denotes the axial vector coupling. The Figure is taken from ref. [p~2j|. 



upper integration limit over dk equal to k in all terms in equations (|36|j37| ) and if we 
set Q 2 in place of W 2 as the upper integration limit in the integral in eq. ([TOD - 

Equations (^,^) were solved in refs. |12[ [L3| assuming the following simple 
parametrisation of the input distributions: 



A P \ U} (x)=N i (l-x) TH (40) 

with i] Uv = r]d v = 3, r/u = i] s = 7 and r\ g = 5. The normalisation constants iVj were 
determined by imposing the Bjorken sum-rule for Au^ — Ad^ and requiring that 
the first moments of all other distributions are the same as those determined from the 
recent QCD analysis JL7|. All distributions Ap\°\x) behave as x° in the limit x — > 
that corresponds to the implicit assumption that the Regge poles which correspond 
to axial vector mesons, which should control the small x behaviour of g\ [|[ || have 
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their intercept equal to 0. It was checked that the parametrisation fl40|) combined 
with equations ( PH| , P^| , |3^| , |3TD gives reasonable description of the recent SMC data on 



gi (x,Q ) and on g\(x,Q ) fl8fl . In Fig.2 we show the nonsinglet part of gi(x : Q 
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Figure 3: The structure function gx(x, Q 2 ) for Q 2 = lOGeV 2 plotted as the function of 
x. Solid line represents this structure function with the double ln 2 (l/x) terms included 
and the dashed line corresponds to g\ obtained from the LO Altarelli-Parisi equations 
The "experimental" points are based on the NLO QCD predictions with the statistical 
errors expected at HERA || . The Figure is taken from ref . [13| . 



for Q 2 = lOGeV 2 in the small x region jL^j. We show predictions based on equations 
([36|,|TI]) and confront them with the expectations which follow from solving the LO 
Altarelli-Parisi evolution equations with the input distributions at Qq = lGeV 2 given 
by equation(s) (f40D. We also show the nonperturbative part of the non-singlet distri- 



bution g 1 



NS{0) l 



x 



5^/6(1 — x) 3 , where gA is the axial vector coupling. In Fig. 3 we 



show gi(x, Q 2 ) for Q 2 = 10 where we again confront predictions based on equations 
([JB|,[n],|n]) with those based on the LO Altarelli-Parisi evolution equations. We also 
show in this Figure the "experimental" points which were obtained from the extrap- 
olations based on the NLO QCD analysis together with estimated statistical errors 
of possible polarised HERA measurements 0. We see that the structure function 
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gi(x,Q 2 ) which contains effects of the double ln 2 (l/x) resummation begins to differ 
from that calculated within the LO Altarelli Parisi equations already for x ~ 10~ 3 . It 
is however comparable to the structure function obtained from the NLO analysis for 
x > 10~ 4 which is indicated by the "experimental" points. This is presumably par- 
tially an artifact of the difference in the input distributions but it also reflects the fact 
that the NLO approximation contains the first two terms of the double ln 2 (l/x) re- 
summation in the corresponding splitting and coefficient functions. It can also be seen 
from Fig. 3 that the (complete) double ln 2 (l/x) resummation generates the structure 
function which is significantly steeper than that obtained from the NLO QCD analysis 
and the difference between those two extrapolations becomes significant for x < 10~ 4 . 
In Fig. 4 we show the spin dependent gluon distribution which contains effects of 




Figure 4: The spin dependent gluon distribution AG(x, Q 2 ) for Q 2 = lOGeV 2 plotted 
as the function of x. Solid line represents AG(x, Q 2 ) with the double ln 2 {l/x) terms 
included and the dashed line corresponds to the AG(x, Q 2 ) obtained from the LO 
Altarelli- Parisi equations. The Figure is taken from ref. |13| . 

the double ln 2 (l/x) resummation and confront it with that which was obtained from 
the LO Altarelli-Parisi equations. It can be seen that the former exhibits characteris- 
tic x~ Xs behaviour with Xg ~ 1. Similar behaviour is also exhibited by the structure 
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function g{(x, Q 2 ) itself. 



To sum up we have presented theoretical expectations for the low x behaviour of 
the spin dependent structure function g%{x, Q 2 ) which follows from the resummation of 
the double ln 2 {l/x) terms. We have also presented results of the analysis of the "uni- 
fied" equations which contain the LO Altarelli Parisi evolution and the double ln 2 (l/x) 
effects at low x. As the first approximation we considered those double ln 2 (l/x) effects 
which are generated by ladder diagrams. The double logarithmic effects were found to 
be very important and they should in principle be visible in possible HERA measure- 
ments (cf. Fig. 3). 
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